Fermionic determinant for dyons and instantons with nontrivial holonomy 
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We calculate exactly the functional determinant for fermions in fundamental representation of 
5*17(2) in the background of periodic instanton with non-trivial value of the Polyakov line at spatial 
infinity. The determinant depends on the value of the holonomy v, the temperature, and the 
parameter n%, which at large values can be treated as separation between the Bogomolny-Prasad- 
Sommerfeld monopoles (or dyons) which constitute the periodic instanton. We find a compact 
expression for small and large r_2 and compute the determinant numerically for arbitrary r\2 and 
v. 



PACS numbers: 11.15.-q,11.10.Wx,11.15.Tk 



in 
o 
o 

(N 

< 

(N 
> 

O 

o 
in 
o 



CD . 

43 



X 



I. INTRODUCTION 

It is well-known that QCD with light fermions has a chiral symmetry restoration phase transition at the temperature 
T ~ 170 MeV. This phase transition has been studied on the lattice and with the help of various QCD-motivated 
models (see, for example [1, 2], [3]), however, not much is known about its driving mechanism. 

When the temperature is high the effective coupling constant becomes small and the theory can be studied semi- 
classically. From the other side, at zero temperature it is very likely that the instanton - anti-instanton ensemble 
(the instanton liquid) is responsible for the spontaneous chiral symmetry breaking [4-7]. Therefore, in order to get an 
insight into the mechanism of chiral symmetry breaking and restoration at non-zero temperatures, it is natural to con- 
sider an ensemble of the (anti) self-dual fields at finite temperature. A generalization of the usual Belavin-Polyakov- 
Schwartz-Tyupkin (BPST) instantons [8] for arbitrary temperatures is the Kraan-van Baal-Lee-Lu (KvBLL) caloron 
with non-trivial holonomy [9, 10]. These configuration were extensively studied on the lattice in [11-13] for SU(2) 
and SU(3) gauge groups. To construct an ensemble of these configurations theoretically, one should first compute the 
quantum weight of the single KvBLL caloron, or the probability with which this configuration occurs in the partition 
function of the theory. In Ref. [14] this problem was solved for the pure Yang-Mills theory. To extend that result to 
the theory with fermions one has to calculate the fermionic determinant in the background of the KvBLL caloron. 
This calculation is the aim of the present paper. 

Speaking of finite temperature one implies that the Euclidean space-time is compactified in the 'time' direction 
whose inverse circumference is the temperature T, with the usual periodic boundary conditions for boson fields and 
anti-periodic conditions for the fermion fields. In particular, it means that the gauge field is periodic in time, and 
the theory is no longer invariant under arbitrary gauge transformations, but only under gauge transformations that 
are periodical in time. As the space topology becomes nontrivial the number of gauge invariants increases. The new 
invariant is the holonomy or the eigenvalues of the Polyakov line that winds along the compact 'time' direction [15] 
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This invariant together with the topological charge and the magnetic charge can be used for the classification of the 
field configurations [16] its zero vacuum average is one of the common criteria of confinement. 

The general expression for the self-dual electrically neutral configuration with topological charge 1 and arbitrary 
holonomy was constructed a few years ago by Kraan and van Baal [9] and Lee and Lu [10]; it has been named the 
KvBLL caloron. In the limiting case, when the KvBLL caloron is characterized by trivial holonomy (meaning that 
the Polyakov line (1) assumes values belonging to the group center Z(N) for the SU(N) gauge group), it is reduced to 
the periodic Harrington-Shepard [17] caloron known before. In this limit the quantum weights were studied in detail 
by Gross, Pisarski and Yaffe [16]. 



'Electronic address: nik_gromov9mail.ru 
^Electronic address: ssliz91ist.ru 



2 




FIG. 1: The action density of the KvBLL caloron as function of z,t at fixed x — y — 0, with the asymptotic value of A4 at 
spatial infinity v = 0.9nT, v = l.l7rT. It is periodic in t direction. At large dyon separation the density becomes static (left, 
ri2 = 1.5/T). As the separation decreases the action density becomes more like a Ad lump (right, ri2 = 0.6/T). In both plots 
the L,M dyons are centered at zl = —vr\2/2nT, zm = v r\2/2nT, il,m = 2/l,m = 0. The axes are in units of temperature T. 



Besides the neutral self-dual configurations there are charged configurations. The "elementary" soliton with unit 
electric and magnetic charges is the dyon, also called the Bogomolnyi-Prasad-Sommerfeld (BPS) monopole [18, 19]. 
It is a self-dual solution of the Yang-Mills equations of motion with static (i.e. time-independent) action density, 
which carries both the magnetic and electric fields of the Coulomb type at infinity, decaying as 1/r 2 . In the SU(2) 
gauge theory there are in fact two types of self-dual dyons [20] M and L with (electric, magnetic) charges (+, +) and 
(— , — ), and two types of anti-self-dual dyons M and L with charges (+, — ) and (— , +), respectively. Formally, L and 
M dyons are related by a non-periodical gauge transformation. 

The KvBLL caloron of the 5/7(2) gauge group (to which we restrict ourselves in this paper) is "made of one L 
and one M dyon, with total zero electric and magnetic charges. It means that for large 7*12 the KvBLL caloron field 
becomes a sum of L and M dyons fields (r\2 is a parameter of the KvBLL caloron field, having a natural meaning of 
the dyons' separation; it is associated with the usual instanton size by ryi = irp 2 T). One can consider the KvBLL 
caloron as a two-monopole solution. Although the action density of the isolated L and M dyons does not depend 
on time, their combination in the KvBLL solution is generally non-static. When the distance r\2 becomes small 
compared to 1/T the KvBLL caloron in its core domain ^Jx^x^ ~ p reduces to the usual BPST instanton (for explicit 
formulas see [9, 22]). The holonomy remains nontrivial as outside the small core domain A4 tends to a constant value. 

The gluonic quantum weight computed in Ref. [14] in the limit when the separation between dyons r±2 is much 
larger than their core sizes 1/v and 1/v (they are constrained by v + v = 2ttT) has the form: 



'8ir 2 Y {Ae~> E \ 3 ( v \£k ( v 

V^rT \ 

exp[-27rri 2 P"(v)] exp -V {3) P(v) , (2) 



where P(v) = ^ 2 ^i T and the overall factor C is a combination of universal constants; numerically C = 1.0314. A 
is the scale parameter in the Pauli-Villars regularization scheme; the factor g~ 8 is not renormalized at the one-loop 
level. 

N f 

To account for fermions we have to multiply the partition function by FJ Det(iy + irrij), where y is the spin-1/2 

i=i 

isospin-1/2 covariant derivative in background considered, and Nf is the number of light flavors. We consider only 
the case of massless fermions here. The operator iy has zero modes [23] therefore a meaningful object is Det'(iy ) 
— a normalized and regularized product of non-zero modes. In the self-dual background it is equal to (Det(— V 2 )) , 
where V is the spin-0 isospin-1/2 covariant derivative [24]. 

In this paper we calculate logDet(— V 2 ) exactly, find a compact analytical expression for the large r%2 asymptotic, 
find 1/r 12 corrections up to the 5 th order, and the small r%2 asymptotic. These results are analytical and give an 
almost exact answer for all ryi except ryi ~ 1/T. In section X we present our results of numerical evaluation of the 
determinant for arbitrary n% and v, which are consistent with both asymptotics. 

In the forthcoming publication [21] we are going to generalize our results to arbitrary SU(N). The first step in 
this direction has been made in Ref. [22] where the gluonic Jacobian over zero-modes, or the volume form on the 
moduli space, has been computed for the general SU(N) KvBLL caloron (the metrics of the caloron moduli space was 
found earlier by Thomas Kraan [25]). We believe that our results will be useful to construct the ensemble of KvBLL 



3 



calorons at any temperatures and to obtain a better understanding of the chiral symmetry restoration mechanism at 
the phase transition temperature. 



II. THE KVBLL CALORON SOLUTION 



In this section we remind some basic facts about the caloron with nontrivial holonomy just to establish the notations, 
see Refs. [9, 10, 14] for a detailed discussion. We use here the gauge convention and the formalism of Kraan and van 
Baal (KvB) [9], the notations are taken from [14]. 

The key quantity characterizing the KvBLL solution for the SU(2) gauge group to which we restrict ourselves 
in this paper, is the holonomy trL, cq.(l). In the gauge where A 4 is static and diagonal at spatial infinity, i.e. 
^Ij-.oc = * v ^j it is this asymptotic value v which characterizes the caloron solution in the first place. We shall 
also use the complementary quantity v = 2irT — v. Their relation to parameters lu,uj introduced by KvB [9] is 
lu = j^p, u) — = 5 — w - Both v and v vary from to 2nT. At v = 0, 2ttT the holonomy is said to be 'trivial', 
and the KvBLL caloron reduces to that of Harrington and Shepard [17]. Note that the fields in the fundamental 
representation feel the sign of trL, therefore the cases v = and v = 2irT differ when we account for fcrmions. 

We shall parameterize the solution in terms of the coordinates of the dyons' 'centers' (we call the constituent dyons 
L and M according to the classification in Ref. [26]): 

T , - 2UJT\2 

L dyon : Z\ = — - 



M dyon : z 2 



T 
2ujt\ 2 



T 1 

dyon separation : z 2 — z 1 = ri2, |r 12 | = irT p 2 , 

where p is the parameter used by Kraan and van Baal; it becomes the size of the caloron at v — > 0. We introduce the 
distances from the 'observation point' x to the dyon centers, 

r = x — z\ = x + 2ujfi2, r = |f|, 

s = x — z 2 = x — 2ujri2, s = \s\. (3) 

Henceforth we measure all dimensional quantities in units of temperature for brevity and restore T explicitly only in 
the final results. 

The KvBLL caloron field in the fundamental representation is [9] (we choose the separation between dyons to be 
in the third spatial direction, r 12 = r 12 e 3 ): 

< vB = \nl„T S d v log* + l -3> Re [{f,l v - i^ v )(ri + ir 2 )d vX ] (4) 

where Tj are Pauli matrices, ff are 't Hooft's symbols [27] with fjfj — e ai j and f\% v = —fj^ 4 = 5 alJ . "Re" means 
hermitian part of the matrix, and the functions used are 

x = Hi ( e -**Bix< — + e 4 ™^ —) $ = t 

i> \ s r J ' ^ ' 

2 

ip = — cos(27rx 4 ) + chch H shsh, tb = tb-\ — — shsh H shch H shch. (5) 

2rs rs s r 

We have introduced the short-hand notations for hyperbolic functions: 

sh = sinh(sv), ch = cosh(sv), sh = sinh(rv), ch = cosh(rv) . (6) 

One can note that the field (4) is not periodical. However, one can formally make a (time dependent) gauge trans- 
formation so that the resulting field is periodical. It turns out that there arc two incquivalent possibilities to make a 
time dependent gauge transformation: g P = e - 2mx i^ T 3 an( j g A = e 2mx 4 uT 3 ^ Correspondingly, we have two self-dual 
periodical fields 

AF(w,x,x 4 ) = Seniors + ^^log* + ^ Re - iffrjfa + ir 2 ) e -^ »d vX ] (7) 
A${u, x, x 4 ) = -Senium + l -nl u Tzd v log* + U Rc [(fj 1 ^ - i^„)(ri + ir 2 )e^^d vX \ (8) 
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Note that they are related by an anti-periodical gauge transformation gpg\ — e^~" lx ^ 3 . In general, an anti-periodical 
gauge transformation may affect different quantities. For example, the determinant of the ghost operator in the 
fundamental representation Det(— V 2 ) is not invariant under this transformation, whereas the determinant in the 
adjoint representation Det(— D 2 ) is invariant. We have to calculate Dct(— V 2 ) for both fields (7) and (8). Fortunately, 
in the case of one caloron we can limit ourselves to one of the fields by the observation that these fields differ only by 
the exchange of constituent dyons and by the substitution uj —* u>. More precisely, 

Ap(u,Xi,X2,x 3 ,X4,) = s M „ e~™^A^(u>,Xi, -x 2 ,2[u> - ui]r 12 - x 3 ,x 4 )e™^ (9) 

where s is a diagonal matrix such that sn = S44 = +1 and S22 = S33 = — 1. We make a rotation around e±, which 
exchanges the two dyons, and a global gauge transformation. Obviously space rotation and gauge transformation of 
a background do not change the determinant and we can conclude 

Det(-V 2 L4»]) = Det(-V 2 L4>)]) (10) 

This allows us to consider only in what follows. We shall systematically drop the P subscript. 

The first term in (7) corresponds to a constant A 4 component at spatial infinity (A 4 «iv^) and gives rise to the 
non-trivial holonomy. 

In the situation when the separation between dyons r\ 2 is large compared to both their core sizes \ (M) and = (L), 
the caloron field can be approximated by the sum of individual BPS dyons, see Figs. 1 (left). We give below the field 
inside the cores and far away from both cores. 



A. Inside dyon cores 

In the vicinity of the M dyon and far away from the L dyon (rv 3> 1) the field becomes that of the M dyon. It is 
static up to the corrections of the order of e ~ 27rwr as can be seen directly from eq.(7) as e -4 ™ 41 '^ become static in 
this domain. We write it in spherical coordinates centered at z 2 : 

iM _ * T 3 ( \ 1\ A M _.sW>Ti -COS0T 2 



Af = ^ v coth(vs) - - , Ag = v 
4 2 V s J 6 2i sinh(vs) 



M M cosj> ti + sin^ r 2 tan(6>/2) 

K =°' A * =V 2 JS inhM +iT3 ^^' 

In the vicinity of the L dyon the field is substantially time dependent. However, this time-dependence can be 
removed by an anti-periodical gauge transformation. It is instructive to write the L dyon field in spherical coordinates 
centered at z\\ 

A L A = ^ ^+27r-vcoth(vr)^ , A^=0, 

L ._— sin(27ra4 — <ft) t\ + cos( 2 7TX4 — </>)t2 

Ae = tV 2sinh(vr) ' 

L _cos(27ra;4 — (f>) t\ + sin(27rx 4 — 0)r 2 tan(0/2) 

A * = lV 2^nhM lT3 ^r— • (12) 

We shall use the fact that the L dyon field can be obtained from the M dyon field by interchanging v with v and by 
making an appropriate anti-periodical gauge transformation. 

We see that in both cases the L,M fields become Abelian at large distances, corresponding to the (electric, magnetic) 
charges (— ,— ) and (+,+), respectively. The corrections to the fields of M and L dyons are hence of the order of 
l/ri2, arising from the presence of the other dyon. 



B. Far away from dyon cores 



Far away from both dyon cores (rv » 1, sv » 1; note that it does not necessarily imply large separations - the 
dyons may even be overlapping) one can neglect both types of exponentially small terms, O (e~ rv ) and O (e~ sv ). 
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With the exponential accuracy the function \ m e q-(5) is zero, and the KvBLL field (4) becomes Abelian [9]: 

•ir ^(v-- -). 




(ri2 - r + s)(r 12 + r - s) 
(r 12 + r + s)(r + s - r 12 ) ' 

In particular, far away from both dyons, A 4 is the Coulomb field of two opposite charges. 

III. THE SCHEME FOR COMPUTING DET( V 2 ) 

As explained in section I, to find the quantum weight of the KvBLL caloron, one needs to calculate the small 
oscillation determinant, Det(— V 2 ), where V M = + and is the caloron field [9] in the fundamental represen- 
tation. We employ the same method as in [14, 28]. Instead of computing the determinant directly, we first evaluate 
its derivative with respect to the holonomy v, and then integrate the derivative using the known determinant at v = 
or v — 2ir [16] as a boundary condition. 

If the background field A^ depends on some parameter T 3 , a general formula for the derivative of the determinant 
with respect to such parameter is 



d logDet(-V 2 L4]) 



-fd 4 xTr(d r A^J^), (15) 



dV 

where is the vacuum current in the external background, determined by the Green function: 

Jf= (8 a J b d x -S a J b d y +A ac d b +A db S a c )g cd (x, y)\ y=x or simply J M = + gV ^. (16) 
Here Q is the Green function or the propagator of spin-0, isospin-1/2 particle in the given background A^ defined by 

-V 2 x G(x,y)=SW(x-y) (17) 

The anti-periodic propagator can be easily obtained from it by a standard procedure: 

+00 

9(x,y)= (-l) n G(x 4 ,x;y 4 + n,y). (18) 

n— — 00 

Eq.(15) can be easily verified by differentiating the identity logDet(— D 2 ) = Trlog(— D 2 ). The background field A^ in 
eq.(15) is taken in the fundamental representation, as is the trace. Hence, if the anti-periodic propagator Q is known, 
eq.(15) becomes a powerful tool for computing quantum determinants. Specifically, we take V = v as the parameter 
for differentiating the determinant, and there is no problem in finding d v A^(r 7 s,ri2,x ,v) for the caloron field (4). 
We assume positions of dyons fixed for convenience, it differs simply by a global translation from d v A^ with the fixed 
center of mass position. 

The Green functions in the self-dual backgrounds are generally known [29, 30] and are built in terms of the Atiah- 
Drinfcld Hitchin-Manin (ADHM) construction [31] for the given self-dual field. A subtlety appearing at nonzero 
temperatures is that the Green function is defined by eq.(17) in the Euclidean IR space, where the topological charge 
is infinite because of the infinite number of repeated stripes in the compactificd time direction, whereas one actually 
needs an explicitly anti-periodic propagator (18). To overcome this nuisance, Nahm [29] suggested to pass on to the 
Fourier transforms of the infinite-range subscripts in the ADHM construction. We review this program in Appendix A 
of [14], first for the single dyon field and then for the KvBLL caloron. This way, we get the finite-dimensional ADHM 
construction both for the dyon and the caloron, with very simple periodicity properties. The isospin-1/2 propagator 
in R 4 is known to be 

In what follows it will be convenient to split it into two parts: 

G(x,y) = g\x,y)+g s (x,y) 
g s (x,y) = G(x,y), g r (x,y) = J2(T^ n G(x 4 ,x;y 4 + n,y). (20) 

n#0 
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The vacuum current (16) will be also split into two parts: "singular" and "regular" , in accordance to which part of 
the periodic propagator (20) is used to calculate it: 

J» = Jl + Jl- (21) 

Although there is no principle difficulty in doing all calculations exactly for the whole caloron moduli space, at 
the stage of spatial integration of the determinant density we loose the capacity of performing analytical calculations 
for the simple reason that the expressions become too long, and so far we have not been able to put them into a 
manageable form in the general case. Therefore, we have to adopt a more subtle attitude. First of all we restrict 
ourselves to the part of the moduli space corresponding to large separations between dyons (r\2 3> 1). Physically, it 
seems to be the most interesting case. Furthermore, at the first stage we take ri2V, r^v 3> 1, meaning that the dyons 
are well separated and do not overlap since the separation is then much bigger than the core sizes, see Fig. 1 (left). In 
this case, the vacuum current (16) becomes that of single dyons inside the spheres of some radius R surrounding 
the dyon centers, such that ~, = -C R -C ri2, and outside these spheres it can be computed analytically with the 
exponential precision, in accordance with subsection II. B, see Fig. 2. Adding up the contributions of the regions near 
two dyons and of the far-away region, we get dDet(— D 2 )/dv for well-separated dyons. Integrating it over v we obtain 
the fermion determinant itself up to a constant and possible l/ri2 terms. Finally, we compute the l/ri2 corrections 
up to the l/rf 2 term in the Det(— D 2 ), which turn out to be quite non-trivial. 

This is already an interesting result by itself, however, we would like to compute the constant, which can be done 
by matching our calculation with that for the trivial caloron at v = and v = 2nT . It means that we have to extend 
the domain of applicability to ri2V = 0(1) (or ri2V = 0(1)) implying overlapping dyons, presented in Fig. 1 (right). 



IV. NEAR EACH OF THE DYONS 



In this section we calculate the right hand side of eq.(15) in the core domains of the L and M dyons. For M dyon 
this task is nearly solved in [14]. In Appendix A the derivation of the vacuum current in the background is given and 
we have to multiply it by an expression for the gauge field (11) and integrate over M-dyon core domain (see Fig. 2). 
The result is 



(1 — svcoth(sv)) (sinh(2sv) — 6sv) (1 — sv coth(sv)) (sinh(2sv) — 2sv) 



967r 2 s 3 sinh 2 (sv) 

and integrating over the ball of radius R surrounding the dyon we obtain 



24ssinh (sv) 



(22) 



V 



7T 

432 



IE 

12tt 



31 

1447T 



log(v/7r) 
12tt 



f R d V P 











4ns 2 ds 



(23) 



where we denote P(v) = r^s (in fact this function is a one-loop perturbative effective potential [16, 32]) and thus 



R : '> 



R 2 



R 



96 



(24) 



7 



The calculation in the L-dyon background is similar but slightly more difficult. The strategy is to reduce this problem 
to the calculation in the M-dyon background, which is simpler. For this purpose we note that the L-dyon gauge field 
is a gauge transformation of the M-dyon gauge field with the parameter v taken to be equal to v = 2ir — v. We should 
note, however, that this transformation is anti-periodical in time and the anti-periodical Green function becomes a 
periodical one under this gauge transformation. It means that we have to the use the periodical vacuum current in 
the M-dyon background computed in the Appendix A and substitute v with v to get the right answer: 



tr[d v A^J l 



L tL-i 

ii J 



(rvcoth(rv) - l)(sinh(2rv) - 2rv) ^ (rv coth(rv) - l) 3 (sinh(2rv) - 6rv) 

12r sinh 2 (rv) 96n 2 r 3 sinh 2 (rv) 

rvcoth(rv) — 1 rv[2 cosh(rv) + cosh(2rv) — 4rvcoth (^) + 7] — sinh(2rv) 



coth(f) 

Integration leads to the following result 



/' 

Jo 



tv[d v A^]d* X = -— - J± 



32nr 2 sinh 2 (rv) 
31 log(v/7r) 



/0 432 12tt 144tt 12tt 

Finally, substituting (23) and (26) into (15) we arrive at the contribution from the cores 



<9 v logDct(-V 2 )| 



log(v/v) 47Ti? 3 

12tt 3 I 2 



T^pm ( v 
2 I 2 



47rr dr. 



(25) 



(26) 



(27) 



CONTRIBUTION FROM THE FAR REGION 



In this region we can drop the exponentially small terms e~ vs and e~ vr in the r.h.s. of eq.(15). It turns out to 
be a great simplification. The calculations are similarly to the calculations for isospin-1 [14]. We are not repeating 
them here because the result is simple and natural. The full vacuum current can be easily expressed through the 

2—2 

perturbative potential energy P(v) = j^j- The only non-zero component is 



J4 = -? p '^Uh^w) 

Making use of eq.(13) we come to the contribution from the far region 

a v logDet(-V 2 )| far = J d v P 



1/1 1 

- v H 

2 V s r 



d 4 x 



(28) 



(29) 



The integral is taken over the whole space with two holes (see Fig. 2). We can use the symmetry between dyons to 
write 



MH T+ H) 



d 4 x 



MI)Ht*0-iH 2 >""- 8 - > +0 (T 



(30) 



where V is the 3D volume of the system. 



VI. COMBINING ALL THREE REGIONS 



Adding (30) to (27) we get the expression for the v derivative of the determinant for large distances r 12 between 
constituents with the l/ri 2 precision: 



(31) 



This equation can be easily integrated over v up to a constant, which in fact can be a function of the separation 
ri2- 



logDet(-V 2 ) 



far 



PI- -- 

2 J 12 



V + P" I 717,12 vlo S v vl ogv 



12tt 



12tt 



+ /(n 2 ) 



(32) 
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This result is valid with 0(exp(— vri.2), cxp(— vri.2), I/V12) precision, and we have separated a constant = P(n), 
which accounts for normalization to the v = determinant. Since in the above calculation of the determinant for 
well-separated dyons we have neglected the Coulomb field of one dyon inside the core region of the other, we expect 
that the unknown function f{r\2) = 0(\jr\2) + ci where ci is the true integration constant. One should consider 
the ri2-derivative of the determinant to sec this explicitly. We perform it in Section VIII. Our next aim will be to 
find the constant ci . The C(l/ri2) corrections will be found later. 



VII. EXPRESSION FOR THE CONSTANT IN THE DETERMINANT 



To find the integration constant, one needs to know the value of the determinant at v = 0, where the KvBLL 
caloron with non-trivial holonomy reduces to the Harrington-Shepard caloron with trivial holonomy, and for which 
the determinant has been computed by Gross, Pisarski and Yaffe (GPY) [16]. Before we match our determinant at 
v^O with that at v = 0, let us recall the GPY result for the isospin-1/2 determinant. 



A. Det(-V 2 ) at v = 

The v = 0, 2n, periodic instanton is traditionally parameterized by the instanton size p. It is known [16, 33] that 
the periodic instanton can be viewed as composed of two BPS monopoles one of which has an infinite size. It becomes 
especially clear in the KvBLL construction [9, 10], where the size of the M (L) dyon becomes infinite as v — > (2tt), see 
section II. Despite one dyon being infi nitely large, one can still continue to parameterize the caloron by the distance 
ri2 between dyon centers, with p = y/r^/ir. Since the determinant (32) is given in terms of r\i we have first of all 
to rewrite the GPY determinant in terms of ri2, too. Actually, GPY have interpolated the determinant in the whole 
range of p (hence ru) but we shall be interested only in the limit ryi ^> 1. In this range the GPY result reads: 

7 2^ , ._7rrL2 1 , , c , m ( 1 



logDct(-V 2 ) - logDct(-V 2 )| v T=0 + r,^- - - logr 12 + ^ + O ^— j (33) 

where r/ = — 1/2 for anti-periodic boundary conditions (i.e. v = 0, and M-dyon is infinitely large), 77 = 1 for periodic 
conditions (i.e. v = 2n, and L-dyon is large). One can verify that the term linear in d in this expression is exactly 
equal to the term linear in d in eq.(32) when v = for anti-periodic and when v = 2n for periodic boundary conditions. 
In Ref. [14] we have managed to find a constant cq analytically: 



-fE 2 TT 4 k>g7T 



9 3 27 
The zero-temperature instanton determinant is [27]: 



logDct(-V 2 )| vT=0 = Ilog M +llog(^)+a(l/2), (34) 

where it is implied that the determinant is regularized by the Pauli-Villars method and p is the Pauli-Villars mass. 
Thus the isospin-1/2 result for the case of trivial holonomy is 

logDet(-V 2 ) = ry^f- + H + a(l/2) - ^ logTr + 1 logp , (35) 

where 

ad/,,. " + a + is:-£W (36) 

72 6 6 tt z 

We have to match our eq.(32) with eq.(35) at v = or v = 27r, but eq.(32) has been derived assuming ri2 ^> ^, = 
and one cannot formally take the limit v — > in that expression without taking simultaneously 7*12 — > 00. In the next 
subsection we will show that one can actually take this limit and get a correct constant. We will rely on exact results 
to show that eq.(32) is valid for arbitrary values of vri2 and vri2 (but large ri2 3> 1). 
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B. Extending the result to arbitrary values of vri2 



Let us take a fixed but large value of the dyon separation r 12 1 such that both eq.(32) and eq.(35) are valid. 
Actually, our aim is to integrate the exact expression for the derivative of the determinant 

8 V logDct(-V 2 ) = J p{x)d i x, p{x) = -Tr [d v A^} , (37) 

from v = 0, where the determinant is given by eq.(35), to some small value of v <C 1 (but such that vr 12 1), where 
eq.(32) becomes valid. We shall parameterize this v as v = k/r\ 2 <C 1 with k ^> 1. The result of the integration over 
v must be equal to the difference between the right hand sides of eqs.(32,35). We want to show that it is 0{k/r\ 2 )', it 
means that there are no large corrections to eq.(32) at small v that can alter the constant term in the determinant. 
Denoting by p our p with subtracted asymptotic terms we have to show that 

j^ 12 dv J d i xp{x) = O (J-^j . (38) 

In this integration we are in the domain 1/v » 1 and n 2 3> 1 and we can simplify the integrand dropping terms 
which are small in this domain. At this point it will be convenient to restore temporarily the temperature dependence. 
With (3 = 1/T our domain of interest is 1/v 3> (3 and r\ 2 S> (3- Therefore we are in the small-/? range and can expand 
p in series with respect to (3: 

P - ^Po + ipi + O(p ) . (39) 

As we shall see in a moment, only the first two terms can be not small in this range and we need to know only them to 
compute c 2 . It is a great simplification because po,i does not contain terms proportional to e~ w since v = 2itT—v — > oo 
at (3 — > 0, and what is left is time independent. Moreover, what is left after we neglect exponentially small terms are 
homogeneous functions of r, s, r i2 , v and we can turn to the dimensionless variables: 

po(r,s,r 12 ,v) = — p ( — , — ,vn 2 J . (40) 



7"i2 \ri2 r\2 



From the exact expressions we see that pi = 0. The reason is that there is no 1/(3 term in eq.(A6). Actually the same 
phenomenon can be seen for the single dyon in eq.(22). We rewrite the l.h.s. of eq.(38) in terms of the new variables: 

J d * xp{x) = r J± Jd 3 xp + O([3) = r jlF a (vr 12 )+O(f3) , (41) 

where x = x/r 12 is dimensionless. We see that it is indeed sufficient to take just the first two terms in the expansion 
(39) at (3 -> 0. 

Let us try to integrate eq.(41) over v from to k/r 12 even without explicit knowledge of F . We get 

k/r 12 k 

J V -^F (vr 12 )dv = ^ J F (vr 12 )d(vr 12 ) (42) 
o o 

We know that the approximate result eq.(32) is correct with at least C(exp(— k), l/(vn 2 )) precision, so at large k the 
correction to eq.(32) must be a constant, plus a term of the order of l/(vri 2 ) = 1/k. Therefore, Fq must be zero. We 
have checked numerically that it is indeed so (actually this function F is the same as in [14]). 
Thus, we have proved eq.(38). We can now compare eq.(32) with eq.(35) and get the constant: 

c h = % + «(l/2) + ^ log 4. + i log M = \ logM + 361 ° g(27r ;i: 39 ' 7r2 - ^ = \ logM + 0.36584648 (43) 
2 lb VI b 6 zlb 7r z b 

It means that eq.(32) is right in whole range < v < 27r with the accuracy l/ri 2 . 
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VIII. l/ria CORRECTIONS 

In the previous section we have seen that Eq.(32) is valid to all orders in l/(vri 2 ), l/(vri 2 ); however, there are 
other l/?~i2 corrections which are not accompanied by the 1/v, 1/v factors: the aim of this subsection is to find them 
using the exact vacuum current. 

To this end, we again consider the case r i2 3> ^, = such that one can split the integration over 3d space into three 
domains shown in Fig. 2. In the far-away domain one can use the same vacuum current as we had in Section V, as 
it has an exponential precision with respect to the distances to both dyons. In the core regions, however, it is now 
insufficient to neglect completely the field of the other dyon, as we did in Section IV looking for the leading order. 
Since we are now after the l/ri 2 corrections, we have to use the exact field and the exact vacuum current of the 
caloron but we can neglect the exponentially small terms of order e~ ri2 . 

Another modification with respect to section IV is that we find it more useful this time to choose r\ 2 as the parameter 
V in cq.(15). We shall compute the l/r 2 2 and higher terms in 9Dct(— D 2 )/dr 12 and then restore the determinant itself 
since the limit of r\ 2 — > oo is already known. Let us define how the KvBLL field depends on m- As is seen from 
eq.(7) the KvBLL field is a function of r, s, v, r 12 only. We define the change in the separation r 12 — ► r 12 + <ir 12 as 
the symmetric displacement of each monopole center by ±dri 2 /2. It corresponds to 



dr r 2 2 + r 2 - s 2 ds r 2 2 + s 2 - r 2 



dru 4ri 2 r dm Ar\ 2 s 

We shall use the definition (44) to compute the derivative of the caloron field (4) with respect to r 12 . 

Let us start from the M-monopole core region. To get the l/ri 2 correction to the determinant we need to compute 
its derivative in the l/rj 2 order and expand correspondingly the caloron field and the vacuum current to this order. 
Wherever the distance r from the far-away L dyon appears in the equations, we replace it by r = (r\ 2 + 1sr\ 2 cos 8 + 
s 2 ) 1 / 2 , where s is the distance from the M-dyon and 9 is the polar angle seen from the M-dyon center. Expanding 
in inverse powers of m we get the coefficients that are functions of s,cos6>. One can easily integrate over 9 as the 
integration measure in spherical coordinates is 2ns 2 ds dcos 8. Then we integrate over s, the distance to the core of the 
M-monopole. After that we have to add contributions from the L and M dyon cores and from the far-away domain. 
It turns out that contributions from dyons differ by terms dependent on R only: 

n , 1 / 23 7 7T l0g(i? 2 Vv/7T 2 )\ 1 fl 1\ 

8 ri2 logDet corc = — (— + -L + — + BV l ' + - + - (45) 



r{ 2 V72tt 6tt 216 12tt / 12r 3 9 7r \ v 



1 /l l\ 60-7r 4 /l 1 \ / 1 



24rf 2 7r Vv 2 v 2 7 2160rf 2 7r \v 3 v 3 / v l2 

(we drop here powers of R as they cancel with the contribution from the far domain) 

Now let us turn to the far-away domain. Recalling eq.(28) we realize that the contribution of this region is 
determined by the potential energy: 

d log Dot (-V 2 ) 



dr 



12 



f d A xd r 

far 



I— 1 


1 1" 




v- - + - 


u 


r s 



(46) 



= ip"(v/2) fd 3 xd ri2 (---) +-i-P IV (v/2) [d s xd, ' ' ! 



K r sj ' 384' v ' ' J " " ri2 \ r .s 
The integration range is the 3d volume with two balls of radius R removed. We use 

a f 1 , 167ri?2 ^ 8n7: ( R\ 2n+1 

U ~ s J = 4tt - ^-5— + g (2n _ 1)(2n+1) (- J (47) 

1 1\ 4 , 3 _ 327rlog(ri 2 /i2) , 2tt(8 - 9^ 2 ) 807ri? 32^i? 2 3047ri? 3 6 

~ _ ~) d X - 12 + 53 + —3 14 + + U \ i / r 12) 

' °J '12 J 'l2 '12 '12 10 '12 



Adding up all three contributions we see that the region separation radius R gets canceled (as it should), and we 
come to: 

^logDet(-V') = Ip»(I) +5 i-(l„ g (^ 2 /^,-a|! + 2 1 + |) (48) 



12r 3 9 7r V v v J 24rf 9 7T V v 2 v 2 / 2160r? 9 7r V v 3 v 3 



r 6 

'12 
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Integrating over r\i gives 



logDet(-V 2 ) 



1 



12n 2 7T 



log(vvr 2 2 /7r 2 ) 



23tt 2 
72 



37 
~6~ 









) 24r 2 2 7r 







(49) 



+ 



1 



72rJ 2 7r \ v 



1 



1 



60 - 7T 4 / 1 



8640rf 2 7r Vv 3 



1 



12-7T 4 



2880rf 2 7r \ v 4 



12 



where c is the integration constant that does not depend on 7*12 . Comparing eq.(49) with eq.(32) at 7*12 — > 00 we 
conclude that 



c= V 



? log\ 
12tt 



flog^ 
12tt 



+ Ci 



(50) 



and ci is given in eq.(43). We note that the leading correction, logr^/ri^, arises from the far-away range and is 
related to the potential energy, similar to the leading r\i term. This series in I/V12 is asymptotic as the coefficients 
are rapidly increasing. 

Interesting, we have revealed that the contribution from the regular part of the current J r has no corrections 1 / r™ 2 
(we verified that up to the 10 th order) and is determined completely by the far asymptote 



(51) 



It seems that Tr {d ri2 A^ JJj) is a full derivative of a simple expression. 



IX. ASYMPTOTE OF SMALL SEPARATION 



In this domain of the moduli space of the KvBLL caloron it is important to realize that KvBLL caloron is a chain 
of the usual BPST instantons equally separated in time direction with a fixed relative gauge orientation between the 
neighbors. 

When the size of instantons is large compared to their separation 1/T, they overlap and cannot be thought of as 
individual pseudoparticles. It is much simpler to consider the KvBL L calo ron as a bound state of two BPS dyons. 
Small separations correspond to the small instanton size p since p = yri2/i. It means that instantons do not overlap, 
and the KvBLL caloron reduces to the usual BPST caloron. 

Formally, one switches to dimensionless variables r — fp, s = sp and then takes the limit p — > 0. In this limit the 
gauge field of the KvBLL caloron becomes that of the ordinary instanton of vanishing size, plus a constant field. It 
is then clear that the determinant becomes that of the ordinary instanton, plus a contribution from a constant field: 

y+ilog/!+i-log(^) +a(l/2) + (D(p) , a(l/2) = 0.145873 (52) 

One can expand in powers of p further, but it requires knowledge of the exact expression for the vacuum current (16). 
Here we give only the final result. The contribution from the regular part of the current J r for small 7*12 is 

-J d^xTr {d ri2 A, j;) = P" Q \ + O(p) (53) 

and from the singular part 

-J d 4 xTr (d ri2 A„ J*) - + (3vv - 2^ 2 )^- + O(p) (54) 

taking into account these contributions we have 

V + a(l/2) + l - log(Mp) - (7tt 2 - 3™ - 3v 2 )^ + 0(p 3 ) , (55) 

Note that the contribution from the regular part of the vacuum current for small values of r\i is exactly the same as 
it has been for large 7*12 in eq.(51). It suggests that eq.(51) is exact. Unfortunately, we have not been able to prove it 
analytically, but we have verified it numerically. We have found out that eq.(51) is correct for all r 12 and v with the 
precision of at least 10~ 3 . 



logDct(-V 2 ) = 



7T 
12 



logDet(-V 2 ) 



PI-)-- 
2 12 



12 




X. NUMERICAL EVALUATION 



One can hardly expect a compact analytical expression for the determinant at arbitrary separations (or instanton 
sizes) and holonomies as even for the case of the trivial holonomy, studied by Gross, Pisarski and Yaffe [16], it is not 
known. In this section we present our results of the numerical evaluation of the determinant for arbitrary separations 
and holonomies. 

We start from eq.(15) and calculate the exact analytical expression for the vacuum current (16) using (19) and the 
exact expression from [14] for the ADHM construction. The resulting expression is too lengthy for printing but can 
be provided by request. We calculate the trace in the integrand in the r.h.s. of eq.(15) with a particular choice of the 
parameter V. 

We have taken V — r\ 2 . Using the axial symmetry we can reduce the number of integrations to 3. We have 
managed to perform the integration numerically with the precision 10~ 4 , despite the complexity of the expression. 
The numerical data are shown in Fig. 3. One can see that it is consistent with all our analytical results. Following 
[16] we denote 



logDet(-V 2 ) = logDet(-V 2 ) 



v,T=0 



A(v,r 12 ) 



P 



n 
12 



V + P" 



■jrr 12 



(56) 



where logDet(— V 2 )| v T _ Q is given in eq.(34). Note that if eq.(51) is exact, only the singular part of the current J s 

contributes to A(v, ri 2 ), thus A(y, 7*12) must be symmetric under the exchange v <-> v. We have verified this symmetry 
numerically. From (49,54) we see that A(y, n 2 ) has the following asymptotics (valid for all v) 



A(y,r 12 ) 



log(27r) vlogv vlogv 



1 

18 



2 

6 



7T 

216 



log(ri 2 /7r) 



6 12tt 12tt 18 6 216 12 

We find that A(v, r\ 2 ) can be fitted by the following expression 



O 



1 

r\2 



(3vv - 27r 2 )ri2 
72^ 



O 



{r% 2 ) (57) 



A(r, r 12 ) 



1 

12 



log(l 



71Tl2 \ 
3 ) 



r 12 a 



216tt(1 +n 2 ) 



0.00302 ri 2 (a + 9vv) 
2.0488 + r 2 



(58) 



where a = 18vlogv + 18vlogv — 216.611. This expression has a maximum absolute error 5 x 10 



XL SUMMARY 



In this paper we have considered finite temperature QCD with light fundamental fermions and calculated the 
fermionic contribution to the 1-loop quantum weight of the instanton with non-trivial value of holonomy (or Polyakov 
line) at spatial infinity. 

Finite temperature theory is formulated in terms of the Euclidean QCD partition function with (anti-)periodic 
boundary conditions in time with the period (3 — l/T, where T is the temperature. 



13 



The instanton with non-trivial holonomy (or the KvBLL caloron [9, 10] for brevity) is the most general self-dual 
solution of the YM equations of motion with a unit topological charge and the above-mentioned periodicity conditions 
imposed. The solution is parameterized by the holonomy v at spatial infinity, the temperature T = 1/(3, the size r\i , 
center position , and color orientation. When r 12 <C fi the solution reduces to the ordinary BPST instanton with size 
P = \Jr\2T j-K in the domain x 2 ~ p 2 which gives the main contribution to the action density. In the opposite limit 
fi2 S> /3 it can be described as a superposition of two (for N c — 2) properly gauge-combed BPS monopoles with the 
separation n 2 between them. 

One can assume that the effective gauge coupling is reasonably small at the deconfincmcnt and chiral symmetry 
restoration temperatures and that therefore the semiclassical ensemble of the KvBLL calorons makes sense and can 
be relevant for describing those phase transitions. It is therefore important to calculate the quantum weight of an 
individual KvBLL caloron as the first step in the study of the ensemble of calorons. The gluonic contribution to the 
quantum weight has been calculated in Ref. [14]. In this paper we have extended the methods and the result of that 
work to the theory with light fcrmions - by computing the fcrmionic determinant over non-zero modes Bet' (if). The 
zero fermion modes [23] should be taken separately when constructing the ensemble, in the spirit of Ref. [5, 6]; it will 
be considered elsewhere. 

The main result of the paper is the determinant in the fundamental representation in the background field of the 
KvBLL caloron (7): 



ilogDet'^y^)) 



logDet(-V 2 )| vT=0 + A(v/T,r 12 T) 



P ( - I - — T 3 
2/ 12 



r - /y '<i Hr' (59) 



where V is a 3c? volume and v = 2nT — v. 

Since the fundamental representation is sensitive to the sign of the asymptotic holonomy e* vTr3 / 2 , tne f er mionic 
determinant is AttT periodic in v, in contrast to the period of 2irT for the gluon and ghost determinants. We have 
also considered the background field (8) for which one has to exchange v with v in eq.(59). 

The numerical values of the quantity A(v, 7*12) are tabulated in Appendix B; they can be fitted by eq.(58) consistent 
with the asymptotic expansion of this quantity at large dyons' separations, 



A(v,r 12 ) = 



+ 



log(27r) vlogv 
log(vvr 



v log v 1 
+ 



12ri 2 7r 
1 



72rf 2 7r 



1 



1 



12tt 

12 A 2 ) - 

60 



18 
23V 
72 

~4 



1 

6 



7T 

216 



37 

+ 2 7+y 



log(ri 2 /7r) 
12 
1 



1 



8640r? 9 7T Vv 3 



1 



24rf 2 7r 

12-7T 4 



(60) 



2880r? 9 7r V v 4 



1 



12 



and with the expansion at small ri 2 (see eq.(55)), 



r i2 (3vv-2^ 2 ) 3/2 



72tt 



+ 0(r 



12 



(61) 



The generalization of these results to the case of larger groups, N c > 2, will be published elsewhere [21]. 
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APPENDIX A: DERIVATION OF THE VACUUM CURRENT IN THE BPS DYON BACKGROUND 

Here we calculate the isospin-1/2 current (16) in the M-dyon background (11). In Section IV it has been explained 
that the calculation in the L-dyon background is equivalent to the calculation for the M-dyon, but with the periodical 
boundary conditions. To calculate the vacuum current, we use the expression for the ADHM quantity \v) for the 
M-dyon from the Appendix A of [14]. It states 

\v) = e ^V^ T V^ T3 ,/ . ™ ~-cxp[zv(ix 4 + st 3 )} (Al) 

y smh(vs) 

We shall use the explicit expression for the Green function in the ADHM background (19). 
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1. Singular part of the dyon current J£ 

The regularized singular part of the current was in fact already computed in Appendix C of [14] . It was named 
there. The result was 



J r s = 0, 

iv (s 2 v 2 csch 2 (sv) + svcoth(sv) — 2) 
487r 2 s 2 sinh(sv) 



4 = A0^2 a 2 ^U^.A (cOs(0)Tl+SnU»T 2 ), 



iv (s 2 v 2 csch 2 (sv) + svcoth(sv) - 2) 

^ = " AO 2 2 ■ x,r \ ^ sm ^ n - cos r 2 , 

A8tt z s z smh(sv) 

= t (sinh 2 (5v) - s 3 v 3 coth(sv)) _^ 

487r 2 s 3 sinh 2 (sv) 

where r are Pauli matrices, s, and 8 are the spherical coordinates centered at the monopole position. Eq.(A2) 
obviously does not depend on the temperature and the type of the boundary conditions as it is a purely zero- 
temperature object. 



2. Regular part of the monopole current 



We are going to calculate the part of the anti-periodical vacuum current corresponding to 

4tt 2 (x- y„) 2 



(CT^^^l-irlgM, y n = y, y n4 = y 4 + n, (A3) 



n#0 

This derivation is similar to that done for isospin-1 in [14]. We repeat it here because it is more simple for the 
isospin-1/2 case owing to the simple form of the isospin-1/2 Green function in the general ADHM background. 
We define 

■>£ = Jf + - ^ + j; 2 = (a; - d^w ■ 

Let us first consider J^ 1 . We have to compute Q r with equal arguments. Substituting (Al) into (A3) one has 

ST(x,x) = Y(-IT I d / v [ cos M2^)l 2 + S inh(2 a vz)r 3 ] emvz 
J Air 2 n 2 smh(sv) 

To compute the sum in this expression we use the summation formula (note that v < 27r): 



(A4) 



e 4Z ™ z 2 \ z \ 1 / a r~ \ 

Era-^-IT + ^' -2 7 r<.<2., (A5) 



in particular 



, 47r 2 n 2 8tt 2 4tt 12 

n^O 



e izn/3 z 2 



It remains now to integrate over 0. The result is 

(3v 2 - 47r 2 )s 2 - 6vcoth(sv)s + 6 



Q T (x,x) 



96tt 2 s 2 



We now turn to the J^ 2 part of the current, where we have to sum over n the derivative of the propagator. First 
of all we consider derivatives of the trace in (A3). One finds for x = y: 
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{&! - 8%){v(z)\v{z)) = ie — — g-^( Tl sin0 - r 2 cos0) , 

sv 

(d$-d*)(v(z)\v(z))=ie mvx . -, r cosOin cos cj> + r 2 sin 0) , 
v ^ sinn(sv) 
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(5| - a|)<«(«)|«(z)) =-2ie mvz — — -[1 2 cosh(2svz) + t 3 sinh(2svz)] , 

smh(sv) 

(af-a?)M«)Kz))=o. (a?) 

The derivative of the denominator of (A3) is zero for x = y except for the derivative with respect to X4, but in this 
case we have the expression of the form of eq.(A4) with n 3 /4 instead of n 2 in the denominator. Now we can sum over 
n. We use the summation formula 

Ee lzn z 3 z\z\ Z . . . 

Next one has to integrate over z. Combining all pieces we obtain: 

j; = 0, (A9) 

iv(3-3svcoth(sv) + s 2 v 2 ) . 

J = ^22.,/ ; (sm0r 2 +COS0T! , 

v 4S>tt z s z smh(sv) 

iv(3-3svcoth(sv) + s 2 v 2 ) , . , , . 

Jg = — 2 2 . , — r sm0ri - COS0T 2 ), 

48n^s 2 smh(sv) 

vcoth(sv)(3 — 3svcoth(sv) + s 2 v 2 ) 1 — svcoth(sv) 

■ h = 48^ + 125 ZT3 ' 

We have used spherical coordinates. For example, the projection of J onto the direction rig = 
(cos 9 cos 4>, cos 6 sin <j), — sin 9) is denoted by Jg . 

The calculations of the regular part of the periodical vacuum current is very similar to that of the antiperiodical 
one. We give here only the result: 

j;+ = o, 

, sv-2tanh(^) 
V = J l + iv —T R — . y sin(/>T2 + cos0Ti ), 
v v 167rsmh(sv)s 

Tr . sv-2tanh(f ) 
J o + *v— - — . , , \ (sin^n - cos^r 2 ), 
Wir smh(sv)s 



Tr+ rr (l - svcoth (f )) (1 - svcoth(rv)) 1 - svcoth(sv) 
h = h 8^ "* + As n 



3. Total monopole current J M 

Summing up the contributions eq.(A2) and eq.(A9) we get for the anti-periodical boundary conditions 

J r = 0, 



ivfl — svcoth(sv)) 2 , 

J * = \o22- u r s m T 2 + cos n , 

4$,tt z s z smh(sv) 



iv(l- svcoth(sv)) 2 

^ = . 2 2 • ^ \ (sm(0)n -cos(0)t 2 ), 

487rs^ smh(sv) 



1 — sv cothfsvl (1 — svcoth(sv)) 



3 



APPENDIX B: RESULTS OF THE NUMERICAL EVALUATION 
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In this appendix we give the numerical data used to draw Fig. 3. 
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